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Abstract 

A non-perturbative treatment is developed for the dephasing produced by the shot noise of a one- 
dimensional electron channel. It is applied to two systems: a charge qubit and the electronic Mach- 
Zehnder interferometer, both of them interacting with an adjacent partitioned electronic channel acting 
as a detector. We find that the visibility (interference contrast) can display oscillations as a function of 
detector voltage and interaction time. This is a unique consequence of the non-Gaussian properties of 
the shot noise, and only occurs in the strong coupling regime, when the phase contributed by a single 
electron exceeds tt. The resulting formula reproduces the recent surprising experimental observations 
reported in [I. Neder et al., cond- mat/0610634 , and indicates a general explanation for similar visibility 
oscillations observed earlier in the Mach-Zehnder interferometer at large bias voltage. We explore in 
detail the full pattern of oscillations as a function of coupling strength, voltage and time, which might 
be observable in future experiments. 

1 Introduction 

Decoherence, i.e. the destruction of quantum mechanical interference effects, is a topic whose importance 
ranges from more fundamental questions like the quantum-classical crossover to possible applications of 
quantum coherent phenomena, such as sensitive measurements and quantum information and quantum 
computing. In mesoscopic transport experiments, decoherence (also called dephasing) is responsible for 
the nontrivial temperature- and voltage-dependence of the electrical conductance in disordered samples 
(displaying weak localization and universal conductance fluctuations) and solid-state electron interferometers. 

The most important paradigmatic quantum-dissipative models ( "Caldeira-Leggett" [TJ [2J and "spin- 
boson" [31 H]) and many well-known techniques used for describing decoherence assume the environment 
to be a bath of harmonic oscillators, where the fluctuations obey Gaussian statistics. This assumption is 
correct for some cases (e.g. photons and phonons), and generally represents a very good approximation for 
the combined contribution of many weakly coupled fluctuators, due to the central limit theorem. However, 
ultrasmall structures may couple only to a few fluctuators (spins, charged defects etc.), thus requiring models 
of dephasing by non-Gaussian noise. Such models are becoming very important right now in the context of 
quantum information processing [51 15j 17] 15] 19] 110 ] ITT1 1121 IT5 1 113] . 

Moreover, the quantum measurement process itself is accompanied by unavoidable fluctuations which 
dephase the quantum system [15H161H7] , while dephasing itself can conversely be viewed as a kind of detection 
process jTSJlIl]- Therefore, "controlled dephasing" experiments can be used to study the transition from 
quantum to classical behavior, e.g. by coupling an electron interferometer to a tunable "which path detector" 
[201 123 HI HI [H [25], which produces shot noise by partitioning an electron stream [251 [271 1281 1251 150] . 
In previous mesoscopic controlled dephasing experiments the coupling between detector and interferometer 
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Figure 1: Schematic drawing of the models considered in the text: A detector channel with shot noise coupled 
to (a) an interferometer channel or (b) a charge qubit. (c) The interaction potential u(x) defines the phase 
function w(x), whose height gives the dimensionless coupling constant g, see equation (fT8|) . (d) Sketch of the 
time-evolution of the oscillations in (a x (t)), indicating the visibility v(t) as the magnitude of the oscillation 
envelope. In this schematic example, the visibility itself oscillates - this is impossible in models of Gaussian 
noise but a direct feature of the non-Gaussian nature of the shot noise, in the strong coupling regime g > tt 
(see text and following figures). 

was weak, requiring the passage of many detector electrons in order to determine the path. Under these 
conditions, the phase of the interfering electron fluctuates according to a Gaussian random process. 

Recently, a controlled dephasing experiment was performed [3H 122] using an electronic Mach-Zehndcr 
Interferometer (MZI) [33J 131] coupled to a nearby partitioned edge-channel serving as a detector. Its results 
differed substantially from those of earlier controlled-dephasing experiments. The interference contrast of 
the Aharonov-Bohm oscillations, quantified by the visibility v = (I max — / ro in)/(-^max + Imin), revealed two 
unexpected effects: 

(a) The visibility as a function of the detector transmission probability T changes from the expected 
smooth parabolic suppression oc T(l — T) at low detector voltages to a sharp "V-shape" behaviour at some 
larger voltages. 

(b) The visibility drops to zero at intermediate voltages, then reappears again as V increases, and vanishes 
at even larger voltages, thus displaying oscillations. 

As estimated in [3T], three (or even fewer) detecting electrons suffice to quench the visibility. For this 
reason, one suspects that these effects may be a signature of the strong coupling between interferometer 
and detector. Indeed, that coupling has already been exploited to entangle the interfering electrons with the 
detector electrons, and afterwards recover the phase information by cross-correlating the current fluctuations 
of the MZI and the detector [31], even after it has completely vanished in conductance measurements. The 
dephasing in the MZI system is caused by the detector's shot noise, which is known to obey binomial, 
i.e. non-Gaussian, statistics. Thus, earlier theoretical discussions of dephasing in the electronic Mach- 
Zehndcr interferometer, based on a Gaussian environment J35[ [36] [37] [38] [39] [40], are no longer sufficient 
(see [HJ [42] for a discussion of Luttinger liquid physics in an MZI). At the same time, a nonperturbative 
treatment is required, to capture the non-Gaussian effects. Higher moments of the noise become important, 
and dephasing starts to depend on the full counting statistics, which itself represents a topic attracting 
considerable attention nowadays [JSJ [3JH [Hj . The relation between full counting statistics, detection and 
dephasing has been explored recently by Averin and Sukhorukov [44]. There, the dephasing rate and the 
measurement rate were considered, i.e. the focus was placed on the long-time limit, similarly to other 
calculations of dephasing by non-Gaussian noise [5j [7] [TO] H21 [131 [14]. In contrast, we will emphasize the 
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surprising evolution of the visibility at short to intermediate times. 

The main purpose of this paper then is to present a nonperturbative treatment of a theoretical model 
that explains the new experimental results, and provides quantitative predictions for the behavior of the 
visibility as a function of detector bias and partitioning. Furthermore, we will show how the approximate 
solution for the MZI is directly related to an exact solution for the pure dephasing of a charge qubit by shot 
noise, where the time evolution of the visibility parallels the evolution with detector voltage. In conclusion, 
it will emerge that the novel features observed in [35], and the results derived here, are in fact fundamental 
and generic consequences of dephasing by the non-Gaussian shot noise of a strongly coupled electron system. 
As a side effect, this may indicate a solution to the puzzling observation of visibility oscillations in a MZI 
without adjacent detector channel [53] . 

The paper is organized as follows: We first describe dephasing of a charge qubit, being the simpler 
model that can be solved exactly. After introducing the model in section |2~TI we derive the exact solution 
(|2.2p . We briefly discuss the relation to full counting statistics (12.31) . and provide formulas obtained in the 
well-known Gaussian approximation (|2.4[ 12.5[) for comparison, before presenting and discussing the results 
obtained from a numerical evaluation of the exact expression (section 12 . 6|) . In section [3j we then go on 
to introduce a certain approximation that keeps only the nonequilibrium part of the noise and allows an 
analytical discussion of many features, some of which become particularly transparent in the wave packet 
picture of shot noise (|3.2p . In section [3J we briefly contrast the features of our solution with those of the well- 
known model describing dephasing by classical random telegraph noise. The Mach-Zehnder interferometer 
is then described in section [5l by first solving exactly the problem of a single electron interacting with the 
detector (|5.1[) , and then introducing the Pauli principle (15. 2[) . The results are discussed (|5 . 3[) and compared 
against the experimental data (|5.4[) . Finally, we briefly indicate (|5.5[) a possible solution to the puzzling 
visibility oscillations observed in the MZI without detector channel. 

Our main results are: the exact formula for the time-evolution of the visibility of the charge qubit given 
in equation (| 19[) . and the formula for the effect of the nonequilibrium part of the noise on the visibility 
of qubit (I37p or interferometer (|65|) . Their most important general analytical consequences are derived in 
section I3~T] including a detailed discussion of the visibility oscillations. 



2 Charge qubit subject to non- Gaussian shot noise 

Interferometers may be used as highly sensitive detectors, by coupling them to a quantum system and reading 
out the induced phase shift. Here we focus on a setup like the one that has been realized in [53], where 
a double dot ("charge qubit") has been subject to the shot noise of a partitioned one-dimensional electron 
channel. However, we note that the strong-coupling regime to be discussed below yet remains to be achieved 
in such an experiment. 



2.1 Model 

We consider a charge qubit with two charge states a z = ±1. It is coupled to the density fluctuations of 
non-interacting "detector" fermions 

H = H qb + H int + H det , (1) 

with H qh = !<x z , 

H dct = y^e fc dj.d fc , (2) 

k 

and 

JTint = ^±lv. (3) 

This coupling is of the diagonal form, i.e. it commutes with the qubit Hamiltonian, thereby leading only 
to pure dephasing and not to energy relaxation (the populations of the qubit levels are preserved). The 
derivation of the exact expressions to be analyzed below depends crucially on this type of coupling. The 
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fluctuating quantum noise potential V introduced in ([3]) is related to the density of detector particles in the 
vicinity of the qubit, see figure [TJ 



V = I dxu{x)pdct{x) =^^u k 'kd\,dk (4) 
J k',k 

Here u(x) is the arbitrary interaction potential (whose details in a realistic situation will be determined 
by the screening properties of the environment), f>det{x) = ip^ (x)ip(x) is the detector density, and ip( x ) = 
^ k <t>k{x)dk is the expansion in terms of the single-particle eigenstates of the detector. At this point we 
do not yet specify the nature of the detector, as some of the following formulas are valid in general for any 
non-interacting fermion system. However, ultimately the evaluations will be performed for a one-dimensional 
channel of fermions moving chirally at constant speed, representing our model of a "detector edge channel" . 
This was implemented in the integer Quantum Hall Effect two-dimensional electron gas [23l [3TJ [32] and can 
presumably be realized in other one-dimensional electron systems as well (e.g. electrons moving inside a 
carbon nanotube). 

We are interested in describing the outcome of the following standard type of experiment in quantum 
coherent dynamics: Suppose we prepare the qubit in a superposition state of ||) and ||) at time t = 0, and 
then switch on the interaction with the detector electrons. In effect this can be realized by applying a Rabi 
7r/2 pulse to the qubit that is initially in the state ||). During the following time-evolution, the off-diagonal 
element p-\i(t) will be affected by the coupling to the bath, experiencing decoherence. Its original oscillatory 
time-evolution is multiplied by a factor, that can be written as the overlap D(t) = (xiWIxtW) °f the two 
detector states Xl(t) and XtM that evolve under the action of -ffdet and -ffdet + V, respectively. In this way, 
the relation between decoherence and measurement becomes evident [18) : 

D(t) = / e +^deti e -i(Hdet + t>)A _ (5) 

Note that we have set H = 1. This can also be written as 

D(t) = /fexpH / dt'V{t'))) , (6) 



where V(t') is the fluctuating quantum noise operator in the Heisenberg picture with respect to iJdct, and 
T is the time-ordering symbol. The magnitude of this time-dependent "coherence factor" defines what we 
will call the "visibility" 



v = \D{t)\. (7) 

The visibility (with < v < 1) determines the suppression of the oscillations in any observable that is 
sensitive to the coherence between the two levels, e.g. (<x x (i)) = Repn(t). This is depicted in figure Q] (d) . 



2.2 Time-evolution of the visibility: General expressions 

The average in the coherence factor D{t) displayed in equation © is taken with respect to the unperturbed 
state of the detector electrons, which may refer to a nonequilibrium situation. We will assume that this 
initial state can be described by independently fluctuating occupations d\dk of the single-particle states k. 
This assumption covers all the cases of interest to us, namely the equilibrium noise at arbitrary temperature, 
as well as shot noise produced by transmission of particles through a partially reflecting barrier, leading to 
a nonequilibrium Fermi distribution. 

The average ([S]) can be evaluated in a variety of ways, e.g. using the linked cluster expansion applied to 
a time-ordered exponential. However, here we make use of a convenient formula derived by Klich [45] in the 
context of full counting statistics. Denoting as T(A) = J2k< k Ak' k d\,dk the second-quantized single-particle 
operator built from the transition matrix elements Awki we have |45j (for fermions) 

tr[e r ^e r ( B )] =det[l + eV], (8) 
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where A is the operator acting in the single-particle Hilbert space. In general, this formula allows us to 
obtain the average of the exponential of any single-electron operator with respect to a many-particle density 
matrix that does not contain correlations. Indeed, for a state with independently fluctuating occupations, we 
can write the many-body density matrix in an exponential form that is suitable for application of equation 
®: 

p = U k [n k dld k + {l-n k )(l-S k dk)} (9) 

= n fc (i-n fe ) e ^ d "^ ln ^r, (10) 

where < n k < 1 is the probability of state k being occupied (formally it is necessary to consider the limits 
rik — > and — » 1 if needed). Inserting this expression into (jSJ), and defining the occupation number matrix 
nk'k — Sk' k n k , we are now able to evaluate averages of the form 

trle^'.^-'^l^p] = (n fe (l-n fe ))det[l + e li — ^r] (11) 

1 — n 

= det[l + (e li -l)n]. (12) 

The average (O then can be performed by identifying the product of time-evolution operators as a single 
unitary operator, of the form given here. Thus, we find 

D(t) = det[l + (S(t) - l)n], (13) 
where the finite-time scattering matrix (interaction picture evolution operator) is 

S(t) = e i ' ldct *e~ i ''' ldct+ "-'' (14) 

Here u is the interaction from ((4]), and hdet is the single-particle Hamiltonian of the detector electrons 
that is diagonal in the /c-basis: [hdet\ k i k = e fe^fe'fe ■ I n principle, equation (|13p allows us to evaluate the 
time-evolution of the coherence factor for coupling to an arbitrary noninteracting fermion system. 

In practice, this involves calculating the time-dependent scattering of arbitrary incoming fc-states from 
the coupling potential u{x), i.e. determining the action of the scattering matrix. Note that in the case of 
fully occupied states (n k = 1 for all k), the operator n becomes the identity and the determinant reduces to 
the product of scattering phase factors that can be obtained by diagonalizing the scattering matrix. More 
generally, the contributions to D(t) from states deep inside the Fermi sea always only amount to a phase 
factor, which will drop out when considering the visibility v — \D(t)\. 

In the remainder of this paper, we will focus on the specific, and experimentally relevant, case of a one- 
dimensional channel of fermions moving at constant speed Wdot (i-e. using a linearized dispersion relation). 
We will employ plane wave states inside a normalization volume L and first assume a finite bandwidth 
k G [— k c , +k c ]. At the end of the calculation, we will send L and k c to infinity (see below). 

The equation of motion for a detector single-particle wave function ip(x, t) in the presence of the potential 
u(x) is 



which is solved by 



i(d t + v d etdx)ip(x, t) = u(x)ip(x, t), (15) 



ip(x,t) = exp[-i / dt' u(x - v d ett')}ip(x - v dct t,0). (16) 



This corresponds to the action of exp(— i(ftdct + u)t) on the initial wave function. Applying exp(ihdett) 
afterwards, we end up with the same expression, but with ip(x, 0) on the right-hand-side (rhs). In other 
words, the action of the scattering matrix is to multiply the wave function by a position-dependent phase 
factor: 

'S{t)if\ (x) = e~ iw ^ijj(x), (17) 
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where the phase function w(x) is related to the interaction potential, as seen above: 



w(x)= [ dt'u{x-v dct t'). (18) 
Jo 

The phase function is depicted in figure Q] (c) . Two remarks regarding the finite band-cutoff k c are in order 
at this point: As argued above, states deep inside the Fermi sea only contribute a phase factor to D(t). This 
is the reason we obtain a converging result for the visibility v = \D(t)\ when taking the limit k c — > oo in 
the end, whereas D(t) itself acquires a phase that grows linearly with k c . Moreover, strictly speaking the 
relation (fTT|) only holds for states ip(x) that are not composed of fc-states at the boundaries of the interval 
k e [—fc c , +k c ], since otherwise the multiplication by e~ lwtyX ^ will yield contributions that are cut off as they 
fall outside the range of allowed wavenumbcrs. Nevertheless, for the purpose of calculating the visibility, this 
discrepancy between the operators S(t) and e™^ will not matter, as those states only contribute phases to 
D{t) anyway. Thus, we are indeed allowed to write the visibility as 

v{t) = det[l + (S(t) - l)n] = det[l + (e^W - l)n] . (19) 

This is the central formula that will be the basis for all our discussions below. 

We briefly discuss some general properties of the phase function uu(x) and its Fourier transform. The 
matrix elements of w are given by the Fourier transform Wk'k = \w{q = k! — k) of w(x). Thus, they are 
connected to those of the interaction potential u(x) via 

w(q) = / dxe~ lqx w(x) = u(q), (20) 

J iqvdct 

where u(q) = J dxe~ lqx u{x). 

At times v^ett ~3> cr, the phase function w(x) has the generic form of a box with corners rounded on the 
scale a of the interaction potential, see figure [1] (c). The phase fluctuations are then due to the fluctuations 
of the number of electrons inside the interval of length v^ctt. The most important parameter in this regard 
is the height of w(x) inside the interval. This defines the dimensionless coupling strength g, given by 

w(q = 0) 1 f + °° , 
9 = I = / dyu(y). (21) 

The coupling strength determines the contribution of a single electron to the phase (in a regime where we are 
allowed to treat that single electron simply as a delta peak in the density). We will see that all the results 
can be expressed in terms of the dimensionless quantities g, eVt, Va/v^ct, and the occupation probability T 
of states inside the voltage window (as well as the temperature, Ta/vdet, for finite temperature situations). 



2.3 Relation to full counting statistics 

In the context of full counting statistics (FCS) [13], one is interested in obtaining the entire probability 
distribution of a fluctuating number N of particles, e.g. the number of electrons transmitted through a 
certain wire cross section during a given time interval, or the number of particles contained within a certain 
volume. Usually, it is most convenient to deal with the generating function 

X (\) = J2 P » eiXN - ( 22 ) 

JV 

The decoherence function D(t), and thus the visibility v = \D(t)\, are directly related to a suitably defined 
generating function. In the limit a —> 0, the phase function w(x) becomes a box of height g on the interval 
x G [0,Vdott]- Then J dx w(x)p(x) is gN, where N is the number of electrons within the box. Thus we find 
for the visibility 



v = \x(g)\, (23) 

in terms of the generating function x for the probability distribution of particles N. For a finite range a of 
the interaction potential, we are dealing with a fluctuating quantity that no longer just takes discrete values. 
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We emphasize, however, that our main focus is different from the typical applications of FCS, where one 
is usually interested in the long-time limit and consequently discusses the remaining small deviations from 
purely Gaussian statistics. The long-time behaviour of decoherence by a detecting quantum point contact 
has been discussed in 44], where formulas similar to (|19|) appeared. In contrast, we are interested in the 
visibility oscillations as a most remarkable feature of the behaviour at short to intermediate times. In other 
words, the kinds of setups discussed here in principle offer an experimental way of accessing such short-time 
features of FCS, which are otherwise not detectable. 



2.4 Gaussian approximation 

Before going on to discuss the visibility arising from the exact expression (|19[) , we derive the Gaussian 
approximation to the visibility. We will first do so in a general way and later point out that the same 
result could be obtained starting from equation (|19p . If V were a linear superposition of harmonic oscillator 
coordinates, and these oscillators were in thermal equilibrium, then its noise would be Gaussian (i.e. it would 
correspond to a Gaussian random process in the classical limit). This kind of quantum noise, arising from 
a harmonic oscillator bath, is the one studied most of the time in the field of quantum dissipative systems 
(e.g. in the context of the Caldeira-Leggett model or the spin-boson model). In that case, the following 
expression would be exact. In contrast, our present model in general displays non-Gaussian noise, being due 
to the density fluctuations of a system of discrete charges. Thus, the following formula constitutes what we 
will call the "Gaussian approximation" , against which we will compare the results of our model: 



-DGauss(0 = exp 



i(v)t-^ ( dtx ( dt 2 (f5V(t 1 )5V(t 2 )} 



(24) 



>o Jo 

Here SV = V — ( V ) . If we are only interested in the decay of the visibility, we obtain 



VGa.uss(t) = |-DGauss(i)| = CXp 



(25) 



1 

2 Jo Jo 2 

i.e. the decay only depends on the symmetrized quantum correlator. Introducing the quantum noise spectrum 



(SVSV^ =1 dte luJt (^SV(t)SV(0)^ , (26) 

we find the well-known expression 

lnv Gauss{ t) = - [P(SV5V) 2 -^pl. (27) 



This result is valid for an arbitrary noise correlator. Inserting the relation between V and the density 
fluctuations (f5|), we have 

(<*VW) = 2nJ2\u k 'k\ 2 n k (l - n k ,)8{u - {e h , - e k )) (28) 

U k',k 

for the spectrum. In the following, we specialize to the case of one-dimensional fermions moving at constant 
speed (efe = Vdct^)- Then we find: 

InuGaussW = ~ ^2 \ w k'k{t)\ 2 n k (l - n k >), (29) 

k',k 

where the matrix w corresponds to the potential u(x — Vdett') integrated over the interaction time, see 
equations (fT8|) and ([20]) . We could have arrived at this formula equally well by using (fT9]) to write 



u = exp[Retrln(l + (e-™-l)n)], (30) 

and expanding the exponent to second order in w. Equation (|29p will be used for comparison against the 
full results obtained from (|19p numerically below. 
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2.5 Results for the visibility according to the Gaussian approximation 



We now discuss the results of the Gaussian approximation for certain special cases. For definiteness, here 
and in the following, we will assume an interaction potential u(x) of width a, which we will take to be of 
Gaussian form wherever the precise shape is needed: 

u (x) = W-) 2 f u( q ) = Mo£ -("/ 2 ' 2 . (31) 

V TO" 

Other smoothly decaying functions do not yield results that deviate appreciably in any qualitatively impor- 
tant way. The coupling strength (f2"Tj) then becomes 

9= ■ (32) 

Vdet 

At zero temperature, in equilibrium, the evolution of the visibility is determined by the well-known physics of 
the orthogonality catastrophe, which underlies many important phenomena such as the X-ray edge singularity 
or the Kondo effect [46]: After coupling the two-level system to the fermionic bath, the two states |x|) and 
\Xi) evolve such that their overlap decays as a power-law. The long-time limit of a vanishing overlap is 
produced by the fact that the ground states of a fermion system with and without an arbitrarily weak 
scattering potential are orthogonal. The exponent can be obtained from the coupling strength g. We find, 
from flSS) and pOjl. 



fv dct t\ ^ 



T=V=0f 



(t) = const .12=1 (33) 



in the long-time limit. Only the prefactor depends on the precise shape of the interaction u(x). We note that 
the result diverges for a —> 0. The reason is that a finite l/a is needed as an effective momentum cutoff up 
to which fluctuations of the density in the Fermi sea are taken into account. In any physical realization the 
fluctuations will be finite, since then the density of electrons is finite and there is a physical cutoff besides 

1/(T - 

After applying a finite bias voltage, the occupation inside the voltage window is determined by the 
transmission probability T of a barrier (quantum point contact) through which the stream of electrons has 
been sent: rik = T for v^k £ [0, eV]. Then, equation (|29[) yields two contributions, one of which is the 
equilibrium contribution we have just calculated. As a result, the visibility factorizes into the zero voltage 
contribution and the extra suppression resulting from the second moment of the shot noise: 



T=0,V^0r,\ 
^Gauss \ l ) 



Gauss 

The function F{eVt) is given by 



7 ll.\A = 0,^ eX P 



-T(1-T)(1) 2 W 



(34) 



F(eVt) = i — 2 — ( eVt ^ 1 and eVa / v dct < 1) , •_>- j 

V ' \ neVt (eVt > l,a = 0) V ' 

Here the low-voltage (short-time) quadratic rise is independent of the shape of the interaction potential u: 
Only low frequency (long wavelength) fluctuations of the density are important, and thus only the coupling 
constant g enters, being an integral over u(x), see (|2ip . At large voltages there is, in general, an extra constant 
prefactor in front of F that depends on a and the shape of u. However, in contrast to the equilibrium part 
of the visibility (f3"3")l . the limit a — > is finite, and we have evaluated this limit in the second line of (|35p. 

Finally, it is interesting to note that for the present model the fermionic density can be expressed as 
a sum over normal mode oscillators (plasmons) after bosonization. Thus, in equilibrium, the Gaussian 
approximation is actually exact. However, once the system is driven out of equilibrium by a finite bias voltage 
and displays shot noise, the many-body state is a highly correlated non-Gaussian state, when expressed in 
terms of the plasmons, even though it looks simple with respect to the fermion basis, where the occupations 
of different fc-states fluctuate independently. 
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Figure 2: Time-evolution of the visibility (coherence) v = \D(t)\ of a qubit coupled to quantum noise from 
a non-interacting ID electron channel, at zero temperature, after switching on the interaction at t = 0. The 
curves have been obtained by direct numerical evaluation of (|19|) . Left: Decoherence by equilibrium noise, 
for increasing coupling strength g (top to bottom curve), displaying the power-law decay (|33p expected from 
the physics of the orthogonality catastrophe (inset: log- log plot, with dashed lines indicating the expected 
exponents (g/2n) 2 ). Right: Decoherence by shot-noise (at a finite voltage eVa/v^et = 1)- Beyond g = n, 
the visibility displays a periodic pattern, with zeroes and coherence revivals, as a result of the non-Gaussian 
nature of the noise. Dashed lines indicate the Gaussian approximation. The transmission probability of the 
barrier generating the shot noise equals T = 1/2. 

2.6 Exact numerical results for the visibility and discussion 

In the following, we plot and discuss the results of a direct numerical evaluation of the determinant (|19p 
that yields the exact time-evolution of the visibility of a charge qubit subject to shot noise. We focus on the 
zero temperature case, although the formula also allows us to treat thermal fluctuations which lead to an 
additional suppression of visibility. The most important parameters are the dimensionless coupling constant 
g (f2Tj) . the transmission probability T, and the voltage V applied to the detector channel. We will also note 
whenever the results depend on the width a or shape of the interaction potential u{x). 

In figure [21 we have displayed the time-evolution of the visibility v = \D(t)\ as a function of v^tt/a, for 
different couplings. In equilibrium, the curves derived from the full expression (|19p coincide exactly with 
those obtained from the Gaussian theory (|29[) . as expected. The long-time behaviour is given by the power- 
law decay (|33[) arising from the orthogonality catastrophe. However, at finite voltages, with extra dephasing 
due to shot noise, the Gaussian approximation fails: In general, it tends to overestimate the visibility at 
longer times and larger couplings (dashed lines in figure [21 right). The most prominent non-Gaussian feature 
sets in after the coupling g crosses a threshold that is equal to g — it, as will be explained below: For larger 
g, the visibility displays oscillations, vanishing at certain times (for a barrier with T = 1/2) and showing 
"coherence revivals" in-between these zeroes. The zeroes coincide with phase jumps of it in D(t) (see figure 
2]). We will discuss the locations of these zeroes in more detail below. 

Such a behaviour of the visibility can only be explained by invoking non-Gaussian noise. In every Gaussian 

theory, we can employ (e tv ) — e~( v ^ 2 > 0, which directly excludes the behaviour found here (regardless of 
noise spectrum and coupling strength), even though it is still compatible with a non- monotonous evolution 
of the visibility. The simplest available model of dephasing by non-Gaussian noise will compared with the 
present results in section [4] 

In order to obtain insight into the general structure of the solution, we first of all note that the qualitative 
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Figure 3: Evolution of the visibility (density plot) as a function of coupling constant g (vertical) and time 
eVt (horizontal). Visibility oscillations start beyond g — it. Note the unequal spacing between zeroes: The 
first zero occurs at a time eVt\ = 2ir 2 /g for large j> 1, see (|42|) . The spacing of subsequent zeroes is given 
approximately by d(eVt) — 2?t for the regime of couplings considered here. Further parameters: T = 1/2 
and eVa/vdct = 1- The red dashed lines indicate the expected location of visibility zeroes, according to the 
approximation v'(t) for the noncquilibrium part, (|37p . in the limit a — * 0. 
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Figure 4: Effect of the transmission probability and phase evolution, (a) Visibility v = \D(i)\ as a function of 
time for different transmission probabilities T (where g — 5 and eVa/vdct — 1). (b) Corresponding evolution 
of the phase, i.e. the argument of the complex coherence factor D r (t), where the subscript indicates that 
the phase evolution for T = 1 has been subtracted as a reference. 

features (in particular the zeroes of the visibility) depend only weakly on the width a or shape of the 
interaction potential. In fact, these features are due to the non-equilibrium part of the noise, and the 
Gaussian approximation suggests that the limit a — > is wcll-dchned for that part. This is the reason why, 
in the following, we will plot the time-evolution as a function of eVt (instead of Vdott/cr), which is the relevant 
variable. 

In figure [3l we display the timc-cvolution versus the coupling g. The threshold at g = ir is clearly noted. 
Furthermore, the first zero occurs at a time t\ = 2ir 2 /(eVg) which shrinks with increasing coupling (provided 
g 3> 1, see discussion in the next section and equation (I42[) ) . In contrast, subsequent zeroes have a periodic 
spacing that appears to be roughly independent of g, given approximately by 8{eVt) = 2ir for the small 
values of g plotted here. 

Finally, in figure 01 the effects of the transmission probability T on the evolution of the visibility and 
the full coherence factor D{t) have been plotted, indicating the phase jumps obtained at T = 1/2 whenever 
v = \D(t)\ vanishes. 

All of these features will now be analyzed further by restricting the discussion to the effects of the 
nonequilibrium part of the noise. 

3 Nonequilibrium part of the noise 
3.1 General properties 

Within the Gaussian approximation, we noted that the visibility at finite voltages factorized into one factor 
describing the decay due to equilibrium noise and another part describing the effect of nonequilibrium shot 
noise (|34[) . More precisely, the nonequilibrium part of the visibility can be calculated from (|29[) by simply 
restricting the matrix elements of Wk>k to transitions within the voltage window: k, k' G [0, eV/vdet]- This 
has the physical interpretation that only these transitions contibute to the excess noise in the spectrum 
(SVSV^J of the fluctuating potential. In addition, since the equilibrium noise comes out exact in the 
Gaussian theory, we can state that all the non-Gaussian features are due to the nonequilibrium part. 
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Based on these observations, we now introduce a heuristic approximation to the full non-Gaussian theory, 
which works surprisingly well. We will factorize 



v T=0 > v *°(t) « v T = v =°{t) ■ v'(t), (36) 

where v' is the visibility obtained from the full expression (fH?)) after restricting the matrix elements of w in 
the fashion described above. We will denote the restricted matrix as w' . Note that the restricted matrix 
depends on the voltage, in contrast to w itself. 

Since the occupation probability is constant within the voltage window, nu = T , the matrices n and w' 
now commute. This allows a considerable simplification, yielding a visibility that can be written in terms of 
the eigenvalues ipj of the matrix 



w': 



v'(t) =U J \TZ + Te- l ^\. (37) 

Thus, the dependence on the transmission probability has been separated from the dependence on interaction 
potential, time, and voltage, contained within ipj. The results obtained from the exact formula are compared 
against this approximation in figure[5](a). We observe that all the important qualitative features are retained 
in the approximation. Furthermore, the locations of the zeroes come out quite well, while the amplitude of 
the oscillations is underestimated. 

We will now list some general properties of the matrix w' that determines the visibility according to (|37|) : 
(i) The sum of eigenvalues is 



7 ' ipj = trw' = i-eVt. (38) 



3 

(ii) For a non-negative (non-positive) phase function w(x), the matrix w' is positive (negative) semidefmite: 
We can map any wavefunction \tp) to another state \ip') by setting tp' k = ipk only inside the voltage window, 
and ipL = otherwise. Then 

(tp \w'\ tp) = (tp' \w\ tp') = J dx \tp'(x)\ 2 w(x) > 0, (39) 
for a non-negative function w(x), and analogously for a non-positive function w(x). 

(iii) Following the same argument, we can prove that the largest eigenvalue of w' is bounded by the maximum 
of w(x), if maxw(s) > 0: 

(ip \w'\ ip) = / dx \ip' [x)\ 2 w{x) < {ip'\ip') m&xw(x) < maxw(x). (40) 



Analogously the smallest eigenvalue is bounded from below by the minimum (if maxw(a;) < 0). 

At small voltages (short times) (where eVt <C 1 and eVa/v^et C 1), the matrix elements are constant 
inside the voltage window, w' k , k = j^w(q = 0), yielding only one nonvanishing eigenvalue, given by 



<Pi = -g-eVt. (41) 

As a consequence, at sufficiently large g 3> 1, the first zero in the visibility v'(t) will occur when (pi = n, 
implying 

Assuming now that w(x) is non-negative (as is the case in our example, if g > 0), we can immediately deduce 
the following general consequences from properties (i) to (iii) : All of them taken together imply that the rise 
of the first eigenvalue must saturate below maxw(a;) (which approaches g for times v^ ct t 3> o~). Thus, other 
eigenvalues must start to grow, in order to obey the sum-rule. If (and only if) the coupling constant is large 
enough, this may lead to an infinite series of zeroes in the visibility (see below). Therefore, we are dealing 
with a true strong coupling effect. 

We have not found an analytical way of obtaining tpj at arbitrary parameters. However, all relevant fea- 
tures follow from the foregoing discussion and may be illustrated by numerical evaluation of the eigenvalues. 
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Figure 5: Left: Contribution v' of the nonequilibrium part of the noise (i.e. the shot noise) to the suppression 
of the visibility. The full lines have been obtained from the exact result (TK)1) by dividing by the result at zero 
voltage, v v ^°(t)/v v=0 (t). The dashed lines represent the approximation of restricting matrix elements of w 
to the voltage window, see (j37|) . We have set eVa/vdct = 1 in this panel. Middle: Universal curves for the 
eigenvalues tpj of the restricted matrix w' entering the visibility v'(t) in equation (|37p . plotted as a function 
of eVt in the limit a — > 0. Right: Locations of the zeroes in the visibility v'(t) as a function of coupling g 
(compare figures [3] and [6|) . These curves can be obtained from those on the left by taking (2n + 1)tt / (ipj / g) , 
with n = 0, 1, 2, 3, 4, . . . (from bottom to top). The blue dotted line corresponds to the first curve displaced 
by 27r, indicating the periodicity observed for small couplings g. In both panels, the green dashed line shows 
the short-time behaviour ipi = geVt /2ir. 

We note that the limit a — > is well-defined, and we will assume this limit in the following, in which results 
become independent of the shape of the interaction potential. This limit represents a good approximation as 
soon as the time is sufficiently large: Wdct^ 3> o. In that limit, the eigenvalues have the following functional 
dependence: 



Thus the complete behaviour at all coupling strengths can be inferred by numerically evaluating the eigen- 
values once as a function of eVt. This has been done in figure El 

At T = 1/2, the visibility v'(t) will vanish whenever one of these eigenvalues is equal to (2n + 1)tt, where 
n = 0, 1, 2, . . .. Thus, the locations of the zeroes can be obtained from the equation (2n + 1)tt = tfj, or 
equivalently g = (2n + l)Tr/(Pj 9 . The latter equation has the advantage that the rhs is independent of 
g. It has been used in the right panel of figure [6] These curves have also been inserted into figure [3l for 
comparison against the results from the full theory. In particular, the first zero is reproduced very accurately, 
while there are some quantitative deviations at subsequent zeroes. 

The full pattern of the visibility, as a function of interaction time and coupling strength, can become 
very complex due to the large number of lines of vanishing visibility v' . This is depicted in figure [6l 

The analysis in the next section indicates (and the numerical results displayed in figure [6] confirm) that 
the spacing between the subsequent zeroes in the visibility is no longer determined by g, but rather given 
by 2ir/eV (with deviations at higher g). As we will explain in the next section, this corresponds to one 
additional detector electron passing by the qubit during the interaction time. 

3.2 Wave packet picture 

Following Martin and Landauer [47], we introduce a new basis of states inside the voltage window, whose 
width in k-space is set by Ak = eV/vdct- 



%=9-^ J) (eW). 



(43) 




(44) 
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Figure 6: Nonequilibrium part of the visibility, v', as a function of coupling strength g and time (or voltage) 
eVt. The lines of vanishing v' are those plotted in figure The structure will be distorted for a > 0, and 
the full visibility v will be further suppressed at higher couplings by the dephasing due to equilibrium noise, 

see dam. 
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In real space, these states represent a train of wavepackets, spaced apart by Ax = 2-k j Ak = VdctAi, 
corresponding to time bins of duration At — h/eV = 2n/eV. Taking the limit L — ► oo, we have 

^ n ( X ) = &¥(-™Ax) sinc( ^ (;c _ nAx))i (45) 
V 2TT 2 

where sinc(y) = sin(y)/y. These packets move at constant velocity, ip n (x, t) — ip n (x — Vd e tt). Their advantage 
is that they are localized in space and therefore well suited for calculating matrix elements of the phase 
function w(x) (or its restricted counterpart w'). From 



(VvHV<«>= E e^Sw (46) 

' ' fe,fc'e[o,Afc] 



we find 

Afc 



1*1 1>») = ^ / +1 dqe^ n+n '^G n ,- n (\q\)w(qAk), (47) 



-l 



where 



, ■ (48) 

1 — \q\ lorn = rt 

These formulas enable a very efficient numerical evaluation. In general, for large \n\ and \n'\ the corresponding 
wave packets lie outside the range of w(x) and therefore the corresponding matrix elements are small and 
can be neglected. In the limit of small voltages, the matrix w' is already diagonal in this basis (compare 
HD): 



(lp„> \w\ Ipr, 



gevt 

2tt 



Sn'oSnO- (49) 



Moreover, the wave packet basis permits a very intuitive interpretation of the results in the Gaussian ap- 
proximation: First, let us assume that the number of contributing wave packets is large, N = eVt/2n 1. 
The wave packets are orthonormalized, and the phase function w(x) is smooth on the scale Ax — v^ctAt = 
v2ir/eV of these packets (for sufficiently large voltages). As a result, we find that the matrix w is diagonal 
in this basis, up to terms of order l/N: 

(lp n < | tf> | 1p n ) PS 6 n 'nW(x = nAx) + O Cjfj • ( 50 ) 

Therefore, in any sum over the eigenvalues ipj, these can be approximated by the values of w taken at the 
centers of the wave packets. Assuming further that the coupling g <§; 1 is weak, each of the (fj is small. This 
allows us to expand the visibility reduction due to nonequilibrium noise: 

v'(t) = Ii 3 \n + Te-^\ « 1 - + ■■ ■ (51) 

3 

That is the expected result, which has the form of "phase diffusion" , with a contribution from the variance of 
the phase shift exerted by each detector electron. In the limit of a <§C eVt, we get ifj = g for approximately 
N wave packets, and zero otherwise. Then we reproduce equation (|35j) in the long-time limit: 

Y>-Ay = ±," = fV. (52) 

3 

More generally, for any shape of w(x) we can replace 

~ E & n 1*1 ^™) 2 ~ ^ J dxw 2 (x). (53) 
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Defining the effective width of w{x) as 



1 2 



[ f dx w(xj] ' 

u = u r , ;/\ , (54) 

J dxw z (x) 

we can set the totai number of wave packets to be iV = I c r/ Ax. Defining the average phase shift induced 
by a single detector electron as (p = J dx w(x)/(NAx) = (geVt/2ir)/N, we can write 

£>f *Nf. (55) 
o 

Note that in the case of a constant w(x), we have I c r — Vdctt, and therefore N = eVt and Cp — g, so we are 
back to {55]). 



4 Comparison with dephasing by classical telegraph noise 

In this section, we compare and contrast the results we have obtained in the rest of this paper against the 
simplest possible model displaying non-Gaussian features in dephasing: pure dephasing of a qubit by classical 
random telegraph noise. There, visibility oscillations are observed once the coupling strength becomes large 
as compared to the switching rate of the two-state fluctuator producing the telegraph noise. In the limit of a 
vanishing switching rate, the visibility in that model is the average of two oscillatory phase factors evolving 
at different frequencies, corresponding to the two different energy shifts imparted by the two-state fluctuator: 



^tel.noise^ = |^ _ ^ + ^-iSut^ (56) 

If the occupation probability of the two fluctuator states is p = 1/2, this leads to visibility oscillations 
i; tcLnolsc (i) = |cos((5o;t) | , roughly similar to those found in our full quantum theory of dephasing by shot 
noise. The decaying envelope of these oscillations is then produced by a finite switching probability. 

It is instructive to set up a rough correspondence between that simple model and the one considered 
here, and see how far it takes us (and where it fails): According to the well-known semiclassical picture of 
binomial shot noise [3D], during the time-interval At = h/eV = 2ir/eV a single detector electron arrives 
with a probability T. It imparts a phase shift g within our model. Thus, the fluctuator probability p would 
equal T, the mean time between telegraph noise switching events would be taken as At, and the frequency 
difference 5cu would have to be set equal to g/ At. This analogy partly suggests the right qualitative behaviour, 
namely a threshold in g that is independent of voltage (independent of At) . This threshold turns out to be 
g = 1 in the telegraph noise model, and for larger g the visibility oscillates with a period Aw At/ \J g 2 — 1 (if 
p = T = 1/2). Although this correctly suggests that the first zero occurs at a position t oc l/(eVg), it predicts 
all subsequent zeroes to occur at the same period, which is not compatible with the actual behaviour (see 
figure [3J. These discrepancies are not too surprising, since the two models certainly differ even qualitatively 
in the following sense: In random telegraph noise, the switching occurs in a Markoff process, i.e. without 
memory. In contrast, in the semiclassical model of binomial shot noise the electrons arrive in a stream of 
regularly spaced time-bins of size At = h/eV. We have not found any reasonable way of incorporating this 
fact into a simplified semiclassical model, since it is unclear how to treat 'fractional time-bins' within such 
a model. 



5 Electronic Mach-Zehnder interferometer coupled to a detector 
edge channel 

In this section we will show how to explain the surprising experimental results that have been obtained 
recently in a strongly coupled "which-path detector system" involving a Mach-Zehnder interferometer coupled 
to a "detector" edge channel. We will present a nonperturbative treatment that captures all the essential 
features due to the non-Gaussian nature of the detector shot noise. Our approximate solution for this model 
is directly related to the exact solution of the simpler charge qubit system discussed above. 
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A simplified scheme of the experimental setup is presented in figure Q] (see [3TJ [35] for a detailed expla- 
nation). Both the MZI and the detector were realized utilizing chiral one-dimensional edge-channels in the 
integer Quantum Hall effect regime. The MZI phase was controlled by a modulation gate via the Aharonov- 
Bohm (AB) effect. The additional edge channel was partitioned by a quantum point contact, before traveling 
in close proximity to the upper path of the MZI, serving as a "which path" phase-sensitive detector [23 . For 
a finite bias applied to the detector channel, the Coulomb interaction between both channels caused orbital 
entanglement between the interfering electron and the detecting electrons, thereby decreasing the contrast of 
the AB oscillations. This contrast, quantified in terms of the visibility v — (/ max — -fmin)/(^max + Anin), was 
measured as a function of the DC bias V applied at the detector channel, and of the partitioning probability 
T of the detector channel. 

As noted already in the introduction, two new and peculiar effects, which will be explained here, were 
observed in these experiments: 

• Unexpected dependence on partitioning: The visibility as a function of T changes from the expected 
smooth parabolic suppression oc T(l — T) at low detector voltages to a sharp "V-shape" behaviour at 
some larger voltages, with almost zero visibility at T = 1/2 (see Fig. 3 in [3"2"]). 

• Visibility oscillations: For some values of the detector QPC gate voltage (yielding T ss 1/2), the 
visibility drops to zero at intermediate voltages, then reappears again as V increases, in order to vanish 
at even larger voltages (see Fig. 4 in [32]). For some other gate voltages it decreases monotonically 
(see Fig. 2 in [31]) • 

In [32) we showed that a simplified model involving a single detector electron can provide a qualitative 
explanation for the experimental results listed above. However, it has clear shortcomings, both quantitative 
and in terms of the physical interpretation. The natural reason for these shortcomings is that detection in 
the experiment is due to a varying number of electrons, not just a single one. Then two questions arise: (i) 
How many electrons dephase the MZI as the detector voltage increases, and (ii) how much does each electron 
contribute to dephasing. These questions will be answered by the following model. 

5.1 Solution of the single-particle problem 

The main simplifying assumption in our approach will be that it is possible to treat each given electron 
in the Mach-Zehnder interferometer on its own, as a single particle interacting with the fluctuations of the 
density in the detector channel. Making this assumption is far from being a trivial step, as it effectively 
neglects Pauli blocking, and we will have to comment on it in the next section. For now, however, let us 
define the following model as our starting point: 



Here x and p = —id x are the position and the momentum operator, respectively, of the single interfering 
electron under consideration (traveling in the upper, interacting path of the interferometer). We have 
linearized the dispersion relation, keeping in mind that the interferometer's visibility will be determined 
by the electrons near the MZ Fermi energy, traveling at a speed wmz- The Fermi energy itself has been 
subtracted as an irrelevant energy offset, and likewise the momentum is measured with respect to the Fermi 
momentum. The Aharonov-Bohm phase between the interfering paths would have to be added by hand. 

We thus realize that the situation is analogous to the model treated above, involving pure dephasing of 
a charge qubit. The two states of the qubit correspond to the two paths which the interfering electron can 
take. The following analysis explicitly demonstrates this equivalence and arrives at an expression for the 
visibility which is the analogue of equation ([6]) . The only difference will be the replacement of i»det by the 
relative velocity v^ct — ^mz, which can be understood by going into the frame of reference of the MZ electron. 

Let us now consider the full wave function |* to tai(*)) = e+^^e"^* |*tot«a(0)) of MZI and detector, 
expressed in the interaction picture with respect to -ffdet- One can always decompose the full wave function 
in the form l^total^)) = / dx \x) ® \tp(x,t)). Here, we focus on the projection \ip(x,t)) = (x l^totaiC*)) on to 
the MZI single-particle position basis. This is a state in the detector Hilbert space, with {ip(x, i)| ip(x,t)) 
giving the probability of the MZI electron to be found at position x. It obeys the Schrodinger equation 




(57) 
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-ivMZ-7^ + V(x,t) 



|#M)>, 



(58) 



where the fluctuating potential V(x, t) e 
to -ffdct ■ The exact solution of equation 



J dx' u(x' — x)pdet(x' ,t) is in the interaction picture with respect 
that follows from the Hamiltonian (fST)) reads: 



\ip(x, t)) = Texp 



-i / dt'V(x-v MZ t',t-t') 



\tp(x - v M zt,0)} 



(59) 



Thus, at a given space-time point (x, t), the "quantum phase" in the exponent is an integral over the values 
of the fluctuating potential at all points on the "line of influence" (x',t') with x — x' = v^zit — t'). If the 
potential were classical, the exponential would represent a simple phase factor. Here, however, the interfering 
electron is not only acted upon by the fluctuations but also changes the state of the detector. The state 
\ip(x, t)) contains all the information about the entanglement between the MZI electron and the detector 
electrons. 

We will now determine the visibility resulting from the interaction between interferometer and detector 
channel. At the first beam splitter, the electron's wave packet is decomposed into two parts, one of them 
traveling along the lower (I) arm of the interferometer, the other one traveling along the upper (u) arm. 
These are described by states \ipi(x,t)) and \ip u (x,t)), respectively, which obey the Schrodinger equation 
given above, albeit in general with a different noise potential for each of them. The visibility is determined 
by the overlap between those two states, taken at the position x = VMzt of the second beam splitter (where 
t is the time-of- flight through the interferometer): 



v = \(Mx,t)\Mx,t))\ ■ 



(60) 



Taking into account that \ip^ u (0,0)) = |\&det) is the detector's initial unperturbed state, and realizing that 
the interaction takes place only in the upper arm, we find that the visibility is determined by the probability 
amplitude for the electron to exit the MZI without having changed the state of the detector [T^ : 









/#det 


Texp 








JO 



dt'V{x - Vu. Z t',t- t') 



del 



(61) 



The initial detector state l^dct) itself is produced by partitioning a stream of electrons. 

The last step consists in representing equation (j6Tj) as an expectation value of a unitary operator: 



del 



- (<I> 



del 



(62) 



where $ is defined as the operator in the exponent of (|61|) . We have been allowed to drop the time ordering 
symbol because the density fluctuations in the one-dimensional detector channel are described by free bosons: 
[p(x, t), p(x', t')} is a purely imaginary c-number. The time-ordered exponential is by definition a product 
of many small unitary evolutions sorted by time. Hence, using repeatedly the Baker-Hausdorff formula 
e A e B = e A+B e^- A - B ^ 2 , which holds since [A, B] commutes with A and B in this case, we can collect the 
operators at different times into the same exponent. The remaining c-number exponent only contributes a 
phase, so it does not lead to a reduction in the visibility and we can disregard it. 

The phase operator $ in (|62[) is therefore a weighted integral over the density operator: 



dt'V(x-v M zt',t-t') = 



dx' w(x')p dct {x')dx' = Wk'kdj/dk- 

k,k' 



(63) 



The phase function w(x) is the one that has been introduced before, in equation (|18|) . with the exception 
that the detector velocity has to be replaced by the relative velocity: t>d c t i"dct — "mz- It can be viewed as 
a convolution of the interaction potential u(x) with the "window of influence" of length I = |(t>dct — "Mz)i| 
defined by the traversal time t and the velocities. 



18 



5.2 Approximate treatment of Pauli blocking 

The loss of visibility is due to the trace a particle leaves in the detector 18J. If the detector (or, in general, 
the environment) is in its ground state initially, this means that the detector has to be left in an excited state 
afterwards. Energy conservation implies that the energy has to be supplied by the particle itself. This is no 
problem if the particle starts out in an excited state. An example is provided by a qubit in a superposition 
of ground and excited state, which can decay to its ground state by spontaneous emission of radiation into 
a zero-temperature environment. 

However, in electronic interference experiments such as the one considered here, we are interested in 
the loss of visibility with regard to the interference pattern observed in the linear conductance. At zero 
temperature, this implies we are dealing with electrons right at the Fermi surface which have no phase space 
available for decay into lower-energy states, due to Pauli blocking. Only an environment that is itself in a 
noncquilibrium state (e.g. the voltage-biased detector channel) can then lead to dephasing. This very basic 
physical picture has been confirmed by many different calculations. While it is, in principle, conceivable that 
subtle non-perturbative effects might eventually lead to a break-down of this picture, we are not aware of 
any unambiguous and uncontroversial theoretical derivation of a suppression of linear conductance visibility 
at zero temperature, for an interferometer coupled to an equilibrium quantum bath. 

The main difficulty in dealing with an electronic interferometer coupled to a quantum bath thus lies in 
the necessity of treating the full many-body problem. Any model that considers only a single interfering 
particle subject to the environment will miss the effects of Pauli blocking, and thereby permit unphysical, 
artificial dephasing by spontaneous emission events that would be absent in a full treatment. In [381 140] . 
it was shown how to properly incorporate these effects into an equations-of-motion approach similar to the 
one described above (with the fermion field tp(x,t) taking the role of the single-particle state \ip(x,t)}). The 
main idea was that the state of the detector, and therefore the noise potential V, will itself be influenced 
by the density in the interferometer, leading to "backaction terms" (known from the quantum Langevin 
equation for quantum dissipative systems) that ultimately ensure Pauli blocking. However, in order to be 
able to solve the equations of motion of the environment, it was crucial to assume Gaussian quantum noise, 
and even then the solution for the visibility was carried out only to lowest order in the coupling. Thus, 
this approach is not feasible for the present problem, where we want to keep non-Gaussian effects in a 
fully nonperturbative way. Nevertheless, the underlying intuitive physical picture remains valid: If both the 
interferometer and the detector are near their ground states, the interfering electron will get "dressed" by 
distorting the detector electron density in its vicinity, but this perturbation is undone when it leaves the 
interaction region. Therefore no trace is left and there is no contribution to the dephasing rate. 

We therefore resort to an approximate treatment (applicable to the zero-temperature situation), suggested 
by the general physical picture described above. We will continue to use the single-particle picture for 
the interferometer, but keep only the noncquilibrium part of the noise, thus eliminating the possibility of 
artificial dephasing for the case when the detector is not biased. In fact, within a lowest-order perturbative 
calculation, this scheme gives exactly the right answers: Firstly, dephasing by the quantum equilibrium noise 
of the detector channel is completely eliminated by Pauli blocking, as follows from the analysis of [38"1 140] 
(at T = 0, for the linear conductance). Secondly, the remaining noncquilibrium part of the noise spectrum, 
corresponding to the shot noise, is symmetric in frequency, and thus equivalent to purely classical noise 
whose effects are not diminished by Pauli blocking (see discussions in [101 HH] ). 

The analysis in section lXTl has demonstrated that all the non-Gaussian features are due to the noncquilib- 
rium part which we retain. Therefore, we expect that the present approximation should be able to reproduce 
the novel features observed in the experiment, which is confirmed by comparison with the experimental data. 
We emphasize once more that it is crucial to supplement the single-particle picture by taking care of the 
Pauli principle afterwards. 

Thus, we shall restrict the matrix elements in equation (|63|) to the voltage window, replacing Wk'k by the 
restricted w' k , k , according to the notation introduced in section [5TT1 All that remains to be done to calculate 
the visibility is diagonalizing the operator $, which is achieved by switching to the basis of eigenstates of w' , 

* = $>ic&, (64) 

3 

where ifj are the eigenvalues and Cj is the annihilation operator for eigenstate j of w' . The occupation 
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Figure 7: Dependence of the visibility v' on the partitioning probability T of the detector current, for 
different voltages, using equation (|55|) . The "V-shape" is clearly observed. Inset: Visibility at T = 1/2 as a 
function of eVt. The locations corresponding to the curves in the main plot are indicated. Other parameters: 
g = 5 and temperature T = 0. 

operators CjCj fluctuate independently, and all states j have the same occupation probability T, just like 
the states in the original basis. This is a consequence of the occupation matrix being proportional to the 
identity matrix, as pointed out near equation (|37p . Therefore, equation (|63p reduces to 

= U j \ll + Te- i ^\. (65) 

This formula is our main result for the visibility of the MZI, valid at zero temperature. It gives a closed 
expression for the reduction of the interference contrast in the AB oscillations of the MZI, as a function of 
detector bias and partitioning probability T. It has been calculated nonperturbatively within the approxima- 
tion discussed above, i.e. employing a single-particle picture for the interfering electron and simultaneously 
retaining only the noncquilibrium part of the detector noise. At any given detector voltage V, there exists a 
basis of states in the detector, which, when occupied, contribute to the MZI phase by different amounts ipj. 
These occupations fluctuate due to the partitioning at the detector beam splitter. The visibility then is the 
product of all those influences. 

5.3 Dependence of visibility on detector voltage and detector partitioning 

The visibility for the Mach-Zchnder interferometer subject to the shot noise in the detector channel may 
thus be calculated in the same manner as the visibility for the charge qubit treated above, if the restriction 
to the nonequilibrium part of the noise is taken into account. The main difference is that in the MZI the 
interaction time t is dictated by the setup. However, in the limit a — > 0, the visibility v' only depends on the 
product eVt. Thus the plots above (figures [5] (b,c) and[6|) also depict the dependence of v' on the voltage at 
fixed time t. 



-ii> 



del 



n>- 



del 
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At small bias voltages, only one eigenvalue is nonzero and grows linearly with detector voltage, according 
to equation (|4Tj) : ipi = geVt/2-K = jV. Thus, the visibility is 

v' = \U + Te-' nV \, (66) 

where the proportionality constant 7 may be measured from the voltage-dependent phase shift obtained for 
the non-partitioned case, T = 1. Equation (f6"6"|) represents the influence of "exactly one detecting electron". 

One can obtain (|66[) as an ansatz, by postulating that exactly one detector electron interacts with the 
interfering electron [32] . Here we obtained it naturally as a limiting case of our full expression. It has to 
be emphasized that this result is highly counterintuitive: Naively, one would assume that each detector 
electron induces a constant phase shift that is set by the coupling strength and does not depend on the 
detector voltage. The voltage V should only control the frequency at which detector electrons are injected. 
However, to the extent that we identify each eigenvalue tpj with one detector electron, we have to conclude 
that this naive picture is wrong. Formally, only a single, very extended detector wavepacket of size cx V^ 1 
interacts with the quantum system, i.e. the charge qubit or the interfering electron (see section |3~H]) . As 
the interaction range is fixed and limited, this means that the phase shift (basically the expectation value of 
w(x) in terms of this wave packet) then shrinks with V. The linear dependence of the phase shift on voltage 
thus may be rationalized by taking into account energy conservation: At lower detector voltages V the phase 
space for scattering of detector electrons gets restricted severely, and thus the effective interaction strength 
is diminished. Likewise, the spatial resolution of this which-path detector becomes very poor, as is apparent 
from the large extent of the wave packet: Detection at a high spatial resolution would prepare a localized 
state that contains a lot of energy, more than is available in the detector-interferometer system. 

Regarding the dependence on the transmission probability of the detector channel (figure [7]), we note 
that there are strong deviations from the smooth dependence exp[— CT(1 — T)] expected for any Gaussian 
noise model (where C would depend on V,t,g,... but not on T). These deviations are particularly strong 
near the voltages for which the visibility becomes zero at T = 1/2. Indeed, if only one eigenvalue contributes 
and is equal to n, equation ([66]) yields a "V-shape" of the visibility, v' = |1 — 2T|, as indicated by the dashed 
line in figure El 

5.4 Comparison with experiment 

In this section we briefly discuss the results obtained by fitting the present model to the experimental data. 
This follows our discussion in [32] . where the reader may find the relevant figures. 

At the outset, we note that the visibility in the real experiment is also suppressed by external low- 
frequency fluctuations, beyond the detector-induced dephasing discussed here. They contribute an overall 
voltage-independent factor that has to be introduced as a fitting parameter when comparing against theory. 

First, we consider the approximation l|66p obtained for low voltages, involving only one detector electron. 
Since the constant 7 was measured, this formula does not contain any free parameters, and can be compared 
directly with the experimental data. As shown in Fig. 3 of [32], it fits very well to the data at low V and 
qualitatively reproduces the novel effects mentioned at the beginning of section [5] In particular, it predicts 
the change from a smooth shape to a V-shape in the dependence of visibility on partitioning probability, as 
well as the non-monotonous behaviour with increasing V. However, according to (|66[) these effects should 
occur when 7V = tt, which does not agree with the experimental observations, where the zero in the visibility 
is shifted to a detector bias that is larger than this estimate by about 40%. Hence at this detector bias ([66]) 
fails to quantitatively reproduce the experimental results. The reason of this discrepancy must be the onset 
of the contributions from other detecting electrons. Once other eigenvalues become slightly non-zero, the 
first one is smaller than jV, because of the sum rule (|38p and the non-negativity of the eigenvalues, (1391) . 
This is clearly apparent in figure [5] The visibility then vanishes at larger values of the detector voltage 
V, in agreement with experiment. At even larger voltages, near tpi « 2-7T, the visibility will have again a 
maximum (coherence revival). However, it will be smaller due to the dephasing by the other detecting states 
(other tpj), again in contrast to the simplified formula ([66]) . These two effects have both been seen in the 
experiment (Fig. 4 in [32]). 

Finally, in [32] , we fitted the experimental data by using, for simplicity, a Lorentzian shape as an ansatz 
for the Fourier transform of the phase function: 
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Here V, has the dimensions of a voltage and turns out to be V* ~ 6.2 fiV. Within this fit, the first eigenvalue 
is if i « 0.87T at V = 9.5fiV, where jV — ir. This implies that almost the full phase shift of ir is contributed 
by a single electron, indicating very strong interchannel interaction. 

5.5 Relation to intrinsic visibility oscillations 

While the earliest implementation of the electronic MZI [55] displayed a rather smooth monotonous decay 
of the visibility with rising MZI bias voltage, this is no longer true in a more recent version [34] . There, 
the visibility displayed oscillations, much like the ones observed here, except they occured as a function 
of MZI bias voltage, in the absence of any detector channel. The present analysis may lead to a possible 
explanation for these initially puzzling observations: The intrinsic intra- channel interaction may cause the 
interfering electrons to be dephased by their own (non-Gaussian) shot noise, if the bias is large enough. We 
note that a similar explanation was put forward in a recent preprint of Sukhorukov and Chcianov |42| , who 
considered a model where two counterpropagating edge channels interacted with each other. Though their 
model is therefore different from ours, we have seen that the visibility oscillations are a generic consequence 
of dephasing by non-Gaussian shot noise, and therefore it is hard to distinguish experimentally (at this point) 
between the different models. 

6 Summary and conclusions 

We presented a nonperturbative approach to the dephasing of a quantum system by an adjacent partitioned 
one-dimensional electron channel, serving as a detector. Our treatment gave an exact expression for the 
time-evolution of the visibility of a charge qubit coupled to such a detector. Moreover, within a certain 
simplifying approximation, it can be used to describe a "controlled dephasing" (or "which path") setup 
where a Mach-Zehnder interferometer is coupled to a detector channel. 

The main features of our results are the following: The visibility may display oscillations as a function 
of time or detector voltage, vanishing exactly at certain points and yielding "coherence revivals" in-between 
those points. This behaviour is only observed if the coupling strength crosses a certain voltage-independent 
threshold, corresponding to a phase-shift of g = ir contributed by a single electron. It is impossible to obtain 
that behaviour in any model of dephasing by Gaussian noise, regardless of the assumed noise spectrum. The 
location of the first zero of the visibility (in detector voltage or interaction time) is proportional to 1 / g for 
large couplings g, while the spacing of subsequent zeroes is approximately independent of g and corresponds 
to injecting one additional detector electron during the interaction time. When plotted as a function of 
detector transmission probability, the visibility differs from the smooth dependence on T(l — T) expected 
for any Gaussian model, rather displaying a "V-shape" at certain voltages. 

All of these features have been observed in the recent Mach-Zehnder experiment [35]. Challenges for 
future experiments include more quantitative comparisons against the theory presented here, as well as 
finding ways of tuning the interaction strength g, to switch between the strong and weak coupling regimes. 
In addition, we hope that the strong coupling physics of dephasing by non-Gaussian shot noise will be seen 
in future experiments involving various other kinds of quantum systems as well. 
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